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1 INTRODUCTION 

Active galactic nuclei (AGNs) are powerful objects in the 
Universe, which are thought to be powered by gravitational 
. ' ener Sy through the accretion of surrounding gas onto super- 
k> massive black holes (SMBHs) located at the center of the 
» I ■ accretion systems. Classically, the accretion paradigm has 
1 been applied to extremely powerful objects like quasars and 
radio galaxies for extragalactic objects. The Eddington ratio 
Lboi/ L-Edd of these objects, where Lboi and Z/Edd are bolo- 
metric luminosity and Eddington luminosity, respectively, is 
~ 10% and thereby geometricall y thin, optically thick dis k 
models, so-called standard disks ( Shakura fc Sunvaevl 19731 ). 



ABSTRACT 

The respective contribution of disk and jet components to the total emission in low 
luminosity active galactic nuclei (LLAGNs) is an open question. This paper suggests 
that 7-rays emitted from electrons accelerated in jets could be a direct diagnostic tool 
for a jet component to the total emission. We demonstrate 7-ray flux from jets based 
on a synchrotron self-compton (SSC) model on the assumption that radio and X-rays 
are dominantly produced from jets in the case of a high state of a nearby LLAGN, 
NGC 4278. We also survey parameter space in the model. Observational properties 
of LLAGNs in radio and X-ray bands allow to constrain physical parameters in an 
emission region. The size of the emission region R is limited to 10 16 cm ^ R ^ 10 175 
cm if the observed tight correlation between radio and X-ray emission originates from 
the same jet component. If the beaming factor of the emission region is close to the 
observed parsec scale jet of NGC 4278 and R ~ 10 16 cm, the 7-rays may be detected by 
Cherenkov Telescope Array, and the jet domination can be tested in the near future. 

Key words: galaxies: active - galaxies: nuclei - galaxies: jets — galaxies: individual: 
NGC 4278 - gamma-rays: galaxies 



have been successfully favored. Recently, the idea that all the 
galaxies have SMBHs at their centers has been common place 
|Kormendv fe Richstone1ll995l : iMagorrian et~ai]|l99g| ). Op- 
tical spectroscopic surveys have revealed that a significant 
fraction of nearby galaxies has also active nuclei, but much 
less luminous than the powerful objects (< 10 40 erg -1 in nu- 
clear Ha lumino siy) , called low luminosity AGNs (LLAGNs) 
|Ho et al.lll997T l. LLAGNs can be furthermore spectroscopi- 
cally classified into Low Ionization Nuclear Emission Region 
nuclei (LINERs), Seyfert galaxies, and transition objects. 
Their Eddington ratios are much lower than the canoni- 
cal value (~ 10%) appropriate for standard disks and could 
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reach ~ 10~ 8 in some cases (e.g.. iHol l|2009i ) ). In addition to 
the low radiative efficiency, correlation among radio luminos- 
ity, X-ray luminosity, and the mass of nuclear black holes 
|Merloni et all 1200a iFalcke et all 120041 ) and the weak fea- 
ture of the big blue bump ( Ho 2008|) indicate optically thin, 
radiatively inefficient accretion flows (RIAFs) in LLAGNs. 

Advection dominated accretion flow (ADAF) models 
|Naravan fc Yilll994 [l995l ; lAbramowicz et al.lll995l ), are a 
kind of RIAFs and have been widely discussed to explain the 
spectral energy distribution (SED) of LLAGNs. ADAF mod- 
els successfully explained most parts of the SED of Sagittar- 
ius A* which is a SMBH located at the center of our Galaxy 
llNaravan et all 1 19951: iManmoto et al.|[l997l : iNaravan et all 
1 19981 : lMahadevanlll99Sl )~ 

ADAF models have been applied to nearby LLAGNs. 
In order to investigate the origin of radiation from LLAGN 
core, radio and X-ray bands are often adopted to avoid possi- 
ble contamination from stellar and dust emission. Follow-up 
observations of optically selected LLAGNs in radio bands 
have detecte d emission dominated with a compact core 
morphology (|Ho fc Ulvestadll200ll : [Nagar et al]|2002f ). De- 
spite their l ow luminosi t y, a l a rger fraction of LLAGNs ar e 
radio-loud (iHo fc Pend I2OOII : iTerashima fc Wilson! 120031 ). 



IDoi et al. I (|2005l ) detected a spectral bump in high-frequency 
radio (~ submillimeter) bands which is a spectral fea- 
ture of an ADAF component. However, observations have 
revealed in radio bands (1-10 GHz) that ADAF compo- 
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nents are not enough to reproduce the total radio emission , 
which indicate the contribution of jets (luivestad fc HollioOll ; 
lAnderson et al.ll2004l : fWu fc Cadl2005l ). These observations 
imply the coexistence of both disk and jet components in 
emission from LLAGNs. 

The respective contribution of disk and j et components 
from r adio to X-ray bands is an open problem. iMerloni et al.l 
(|2003l ) found correlation among radio luminosity, X-ray lu- 
minosity, and black hole mass and showed that an ADAF 
model can reproduce observational data better than a radi- 
ation model dominated by jets by com paring the data with 
theoretical models. On the other hand, iFalcke et alj (|2004l ) 
suggested that radiation is dominated by non-thermal emis- 
sion from jets by a similar analysis. The correlation be- 
tween radio and X-ray luminosities has also implied the com- 
mon non-thermal origin from jets for radio-loud L LAGNs 
llBalmaverde fc CapetH 120061 ; IPanessa et all l2007h . Note 
that IMerloni et al.l ^200% ) also mentioned that considering 
cooling processes of electrons can improve the reproducibil- 
ity of the data in the jet model and synchrotron radiation 
from jets can be res ponsible up to X-rays. Furthermore, 
iNemmen et al.l |201Cf ) demonstrated the SED fittings of ob- 
served LLAGNs and showed that observed X-ray data can 
be reproduced by both a disk-dominated and jet-dominated 
radiation models without inconsistency at present. 

This study suggests that 7-rays emitted from elec- 
trons accelerated in jets can be a direct diagnostic tool 
for a jet component in radio to X-ray bands. When X- 
rays are produced by synchrotron radiation of the electrons, 
the synchrotron photons are also upscattered by the elec- 
trons through inverse compton scattering (ICS) into 7-rays. 
This is so-called synchrotron self-compton (SS C) scenario 
IIMaraschi et al. 1 1 19921 ; iBloom fc Marscherl Il996h ." This sce- 
nario is we ll established for the SED modelings of BL Lac o b- 
jects fe.g jAnderhub et all (|2009l ); lAharonian et all (|2009h V 
We focus on NGC 4278 as an example of LLAGNs and 
demonstrate expected 7-ray flux based on a SSC model on 
the assumption that non-thermal radiation from jets is dom- 
inated. We also survey parameter space in the model and 
constrain physical parameters required in jets. 

NGC 42 78 is a LLAGN wi th z = 0.002 and the distance 
of 16.7 Mpc (|Tonrv et al ll200lh for H = 70 km s" 1 Mpc" 1 , 
which is Hubble constant. This source is sometimes catego- 
rized into LINER or radio-loud LLAGNs. Recent observa- 
tions in radio bands have revealed a two-sided relativistic 
parsec scale jet (/3 ~ 0.75) closely aligned to the line of sight 
(2° ^ ^ 4°), where ft and are the velocity in the unit 
of spe ed of light and the v iewing angle of the jet, respec- 
tively (|Giroletti et allfeoosl ). These observables lead to the 
relativistic beaming factor of the parsec scale jet of 2.6. In 
addition, a monthly time scale varia bility by a factor of 3 to 
5 has been observed in X-ray bands (|Younes et alJlioiol ). In 
general, the time scale of flux variability limits the size of 
emission region. 

This paper is laid out as follows. In Section [2] we de- 
scribe a model of SED originating from electrons accelerated 
in relativistic jets. In Section [3] we fit resultant SED to ob- 
served data below X-rays and estimate 7-ray flux. Moreover, 
combining observational results in radio and X-ray bands, 
we constrain the physical parameters of emission region. Fi- 
nally, we make some discussions and summarize this study 
in Section [4] 



2 MODEL 

Electrons are accelerated in a jet via particle acceleration 
mech anisms, i.e., shock acceleration (|Blandford fc Eichler] 
Il987t ) and emit synchrotron radiation by interactions with 
magnetic field. The synchrotron photons are also upscat- 
tered by the electrons through ICS. For simplicity, this elec- 
tron acceleration region is modeled as a magnetized spherical 
blob with the radius of R, which is moving with the Lorentz 
factor of F in a jet. Once we inject accelerated electrons into 
the blob, we discuss the interactions and emission from the 
blob following a SSC model. In this section, we introduce a 
model of relativistic electron spectrum. The other details of 
our SSC model is described in appendix |A1 

Although statistical particle acceleration mechanisms 
produce a power-law spectrum of electrons, the cooling pro- 
cesses of the electrons make a spectral break. Based on this 
picture and the approximation that electrons are continu- 
ously injected to the blob during the order of R/c, electrons 
accelerated in the blob are assumed to have a broken power- 
law spectrum, 



n 7 



1 + 



Tbi 



(7min < 7 < 7max) , (1) 



where no, 7 are the normalization factor of the number den- 
sity of the electrons and the Lorentz factor of electrons, re- 
spectively. 7 m in, 7max, and 7b r are the minimum, the maxi- 
mum of 7, and and the Lorentz factor at a spectral break, 
respectively. We set 7 m in = 1 throughout this paper, si and 
S2 are the spectral indices below and above 7br- For vari- 
able sources, an electron spectrum is expected to be time- 
dependent in reality, but the time profile of the electron 
injection is uncertain at present. A simple model is instan- 
tenious injection approximat i on dis cussed in literture (e.g., 
iDermer fc Schlickeiserl (| 19931 . l2002h ). On the other hand, 
when particle acceleration at shocks in a jet is considered, 
the shocks propagate over a size of the out-flowing plasma 
a nd can deposite energy cont inuously as also mentioned 
in lDermer fc Sc hlickciscr (j 1993T ) . Thus, continuous injection 
approximation is also plausible in internal shock scenarios. 
It should be constrained by observations which injection ap- 
proximation is closer to reality. In this paper, we adopt the 
continuous injection approximation and the electron spec- 
trum described in Eq[T] 

7max can be estimated by comparing the time scale to 
accelerate electrons and the shortest energy-loss time scale. 
The former time scale is 



race =»f-=6x lO 2 0F,27rB-r 1 s, 



(2) 



speed of light, respectively, and 8f,2 = 0f/W 2 , 77 = 7/IO 7 , 
B-i = B/10 -1 G. We take Of to be a constant, which is 
equivalent to assuming the diffusion coefficient for acceler- 
ated particles to be proportional to the Bohm diffusion co- 
efficient. 9f ^ 10 is a fairly conservati ve value, but 6f ~ 1 
is als o possible for relativistic shocks dRachen fc Meszaros! 
1 19981) . 

The accelerated electrons suffer from energy-loss via 
synchrotron radiation and ICS. The former time scale is 



3m e c 



8 x 10 3 77 _1 B-i" 2 s, 



(3) 
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where m e , or, and Ub = B 2 /8n are the electron mass, the 
cross-section of Thomson scattering, and the energy density 
of magnetic field, respectively. The time scale of ICS has a 
form similar to r syn in the Thomson limit, 

3m e c (4) 



4<TT7t/rad 



where U Ta a is the energy density of radiation. The shortest 
energy-loss time scale is determined by comparing Ub and 
(7 ra d- If £/ r ad is larger than Ub, Tmax and 7br are estimated 
by adopting £/ ra d- The flaring of blazars realize [/ ra d > U ayn 
in many cases. However, (7 ra d cannot be estimated before 
the calculation of SED. Returning to a kinetic equation of 
electrons to determine electron distribution, the equation 
is non-linear because ICS energy-loss rate depends on the 
electron distribu tion. This non-linearity changes resultant 
7-ray spectrum JSchlickeiserll2009l ; ISchlickeiser et all l201Ct 
IZacharias fc Schli ckciscr 20l3). Since we do not know [/ ra d 
initially on the spectral assumption of Eq[TJ first of all, we 
calculate SEDs on the assumption of Ub > (7 ra d, which 
means that the shortest time scale is T syn , and then check 
[/rad- Although this assumption will be satisfied in many 
cases, Ub < frad could be realized when R is small. Fortu- 
nately, such cases are not favored by several observational 
aspects as discussed in Section[3] Consequently, 7 max can be 
estimated by r aC c = Tsyn under this assumption as 



7max — 



6ne 
QfotB 



-1/2 



= 4 x W 7 e F2 ~ 1/2 B. 



-1/2 



(5) 



Note that this 7 max always sat isfies the Hil las criterion that 
r g should be smaller than R (|Hillasl H984h in all the cases 
treated in this paper. 

The cooling of electrons via synchrotron radiation 
and/or ICS makes a power-law index be steeper by one, 
S2 = si + 1 above characteristic energy 7b r , assuming that 
accelerated electrons are continuously provided into the 
blob. Since the electrons lose energy during their staying 
in the blob, this characteristic energy can be estimated by 



syn < t csc , where 



3R 



1 X l(fi?16 S. 



(6) 



is the time scale to escape electrons from the blob. Here, we 
assume that electrons escape from the blob by the velocity of 
c/3, which corresponds to the velocity of downstream fluid in 
a shock system under the strong shock limit of a relativistic 
shock. Thus, 



Ybr 



2nm e c 
a T RB 2 



8 x 10 4 i?i 6 _1 -B_ 



We estimate the SEDs of a jet component of NGC 
4278 in the overall energy range based on the SSC model 
under the electron spectrum modeled above. This model 
has five free parameters: Si(s2), B, R, no, and S, where 
S = [r(l — /3cos6)]~ is so-called the beaming factor of the 
blob. For convenience, we change a parameter no into the ra- 
tio of electron energy density in the blob to Ub, V = U b /Ub, 
where 



U e = m e 



, dn e 
d-yy -j— . 



(8) 



R [cm] 


B [G] 


V 


10 14 


3.72 


140.0 


10 15 


0.80 


37.0 


10 16 


0.165 


10.0 


10 17 


0.035 


2.5 



Table 1. Physical parameters adopted in Figs[T]and[2] S = 2.6 is 

assumed in these figures. 



3 EXPECTED 7-RAYS 

In order to estimate 7-ray flux by fitting a spectrum ob- 
ser ved in radio to X-r ay bands, we adopt data processed 
by lYounes et alj (|2010T ). A selection criterion of this data 
set is high angular resolution less than 10 arc-seconds to 
resolve radiation from the core of NGC 4278. The authors 
adopted the data of both Chandra and XMM-Newton ob- 
tained in different periods in X-ray bands and found a 
monthly time scale variability in which flux is changed by 
a factor of 3 to 5. In this paper, we focus on a high state 
observed by XMM-Newton due to following two reasons. 
The first reason is the existence of near-UV data simulta- 
neously observed with the X-rays by an optical/UV moni- 
tor instrument on board XMM-Newton. Simultaneous data 
for variable sources are generally a powerful tool to unveil 
the origin of emission from the sources. The second rea- 
son is that a high state expects larger 7-ray flux, which 
allows us to diagnose the origin of core radiation more eas- 
ily. Since XMM-Newton does not have angular resolution 
as well as Chandra, it could not resolve two X-ray point 
sources close to the core resolved by Chandra. The authors 
estimated X-ray components other than the core as only 
2% by using Chandra data. The data could be well fitted 
by a power-law spectrum without any thermal component 
and confirmed that a spectral index in X-ray bands is not 
contaminated. The photon index defined in the X-ray spec- 
trum in the unit of photons cm" - 2 s" 1 keV -1 is 2.05 ±0.02. 
Assuming that the X-rays are synchrotron radiation from 
high energy electrons above 7b r , the spectral index of the 
electrons is S2 ~ 3.0. The other optical data was processed 
from data observed by Hubble Space Telescope (HST), but 
they are not observed s imultaneousl y . In radio bands, we 
adopt data observed by Nagar et al.l (|200lNGiroletti et al.l 
(|2005h,|Jones et all (119841 ) . which were also used and listed 



Younes et al.l ( 20101 ). The spectral index a ra d defined as 
vh v oc i/~°™ d in radio bands is ~ 0.5, where L„ is lumi- 
nosity per unit frequency, which corresponds to si ~ 2.0 
because si = 2a r ad + 1 for synchrotron radiation. In order 
to produce SED in overall energy range, we fit the observed 
data in radio and X-ray bands by the model. The variability 
of radio emission is assumed to be small, although the radio 
data are not simultaneous data to X-rays. 

The size of an emission blob can be constrained by sev- 
eral discussions. The monthly time scale variability in X-ray 
bands can limit the blob size as 
ScAt obs 



R < 



l + z 



3 x 10 17 <5 3 At obs , 



1 cm, 



(9) 



where £3 = 5/3 and At b Sj _i = At o b s /10 _1 yr is the vari- 
ability time scale. Although ~ 1.5 hours ~ 10% variability 
in X-ray flux is also reported in the term observed by XMM- 
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Figure 1. Spectral energy distributions on the assumptions of 
5 = 2.6 and 8p = 10 2 . The other parameters are listed in Ta- 
ble [T] The sensitivity limi t of Swift BAT 54 months sky sur- 
vey JCusumano^et^l][20l3) is also shown. In addition, the cor- 
responding integral fluxes, F(> E), are sho wn with the integra l 
sensitivity curves of Fermi LAT (5cr. lyrl l|Atwood et al]|2009l ) 
and CTA (5a, lOOhr) ||CTA ConsortiurrfeOlOh in the small panel. 
The cases of R = 10 14 and 10 15 cm are not favored because of 
synchrotron self-absorption (see text). 
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Figure 2. Same as FigQ] but the cases of Op = 10 



Newton, the origin of this variability does not care here, but 
will be discussed afterward. On the other hand, assuming 
that X-rays are emitted from a blob in jets, its size is difficult 
to be smaller than the Schwarzschild radius of the central 
black hole, r sc h. Sin ce the black hole mass of NGC 4278 is 
M BH ~3x 10 s M q (|Wang fc Zhandl2003l ; IChiaberge et all 
I2005D , the blob size is lower-limited as 



R > 7~sch = 



2GM BH 



= 9 x 10"M, 



0,8.5 



(10) 



where M Q ,8. 5 = M B h/10 8 ' 5 M q , G and M Q are the gravi- 
tational constant and the solar mass, respectively. Thus, we 
consider 10 14 cm ^ R ^ 10 17 cm below. 

FigH] shows SEDs calculated on the assumptions of 
5 = 2.6 and dp — 10 2 . The absorption of high-energy 7- 
rays by extragalactic background light is t aken into account 
by using the model of iFranceschini et al.l (|200Sf ). The other 



parameters required for the spectral fits are listed in Ta- 
ble [1] Note that 5 = 2. 6 was derived from a radio im- 
age of a parsec scale jet (|Giroletti et al.l [20051 ). 8 could be 
larger at inner jets (discussed afterward) as radio observa- 
tions have revealed mildly relativistic parsec scale jets in 
blazars co ntrary to the requiremen t of 8 ^ 10 to reproduce 
7-ray flux (|Piner fc Edwardsl|2004 ). The sensitivity limit of 
Swift Burst Alert Tele scope (BAT) 54 months sky survey 
|Cusumano et afl l201fj ) is also shown. In the small panel 
inside the figure, integral fluxes corresponding to the SEDs 
with the integral sensi tivity curve of Fermi Large Area Tele- 
scope (LAT) (5a, lyr) (|Atwood et al.|[2~009l ) and goal integral 
sensitivity of Cherenkov Telescope Array (CTA) (5<r, lOOhr) 
jCTA Consortiurd[2010l ). 

FigH] demonstrates that a smaller blob predicts larger 
7-ray flux. The case of a smaller blob requires larger number 
density of synchrotron photons, i.e., seed photons for ICS, 
to reproduce the observed flux up to X-rays, which leads 
to larger 7-ray flux. To realize the increase of synchrotron 
photons, a strength of magnetic field and/or the number of 
electrons in the blob should increases. However, only the 
number of electrons is essentially a free parameter, since the 
strength of magnetic field is constrained to reproduce the 
observed data as follows. In order to fit the radio and X-ray 
data, a spectral break of synchrotron radiation should be at 

~ 10 15 Hz. Since the typical frequency of synchrotron 
photons generated from electrons with the energy of 7 is 

v c = 0.29-^^^ = 4 x W 5 S 3 R 16 - 2 B-r 3 l 2 Hz, (11) 
1 + 2 47rm e c 

the requirement of ~ 10 lj Hz means 



7ii 6 2 B_i 3 ~ 7 x 10 _1 <5, 



(12) 

from Eqs[7] and 1111 The values of magnetic field listed in 
Table [1] roughly consistent with this estimation. In Table 
[T] rj is larger for a smaller blob, which means that more 
electrons are required for a smaller blob with respect to the 
increase of magnetic field. Also, the electron energy density 
is dominated compared to the energy density of magnetic 
field for S — 2.6. This is a similar situation to the emission 
region of blazars. 

However, the small size of the blob has several prob- 
lems. Firstly, radiation energy density is dominated com- 
pared to the energy density of magnetic field in the blob 
in the cases of R = 10 14 and R = 10 15 cm. In these cases, 
7max and 7b r are determined by ICS, which are unlikely to 
Eqs[5]and[7] However, even including the effects of ICS, the 
model cannot reproduce the observed data self-consistently 
on the assumption of 8 — 2.6. C/ ra d much larger than Ub 
makes the spectral break frequency of synchrotron radia- 
tion be lower. In order to keep i^br to be ~ 10 Hz, U Ta a 
should be reduced by decreasing n e or B. Both of these 
decreases the flux of synchrotron radiation. In the former 
case, B is required to be larger to compensate the decrease 
of synchrotron radiation, but this leads to a break frequency 
lower than that calculated in Fig[T] Although rj needs to be 
larger to reproduce the observed data in the latter case, the 
total I7 ra d is not reduced because ICS flux is proportional 
to n e 2 . The second reason is another spectral break due to 
synchrotron self-absorption. A smaller R produces the spec- 
tral break at a higher frequency, which is ~ 10 11 Hz and 
~ 10 10 Hz for R = 10 14 and R = 10 15 cm, respectively. In 
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R [cm] 


B [G] 


V 


10 14 


5.86 


0.25 


10 15 


1.26 


0.06 


10 16 


0.272 


0.014 


10 17 


0.058 


0.003 



Table 2. Physical parameters adopted in Fig|3]under <5 = 10. 



these cases, radio emission from the blob is not dominated 
at ~ 5 GHz. For well-observed blazars, radio emission at 5 
GHz is already optically thick and therefore radio emissio n 
from the blob is not dom inated fe.g.. lAnderhub et al. (2009); 
lAh aronian et al.l |2009l )). It is not clear at present whether 
this feature is common in radio-loud LLAGNs. However, the 
tight correlation between radio luminosity at 5 GHz and X- 
ray lumin osity implies that th e origin of these radiations is 
the same (|Panessa et al.l [20071 ). Based on this implication, 
R ^ 10 15 - 5 cm is not favored for S = 2.6. 

The integral flux for 7? = 10 16 and 10 17 cm are lower 
than the integral sensitivity curve of Fermi LAT. This is 
consistent with the fact t hat NGC 4278 is n ot listed in Fermi 
LAT lyr source catalog |Abdo et al.ll2010l ). Unfortunately, 
even 10 years observations do not reach the prediction for 
R — 10 16 cm, assuming that the sensitivity curve of Fermi 
LAT is roughly scaled inversely proportional to the square 
root of the total observation time. The integral flux is also 
compared with the goal integral sensitivity of CTA. The 
integral flux for R = 10 16 cm is above the CTA sensitivity 
at ~ 100 GeV. Thus, CTA can test the jet domination of 
the total emission in this source in the case of R ~ 10 16 cm. 

The simultaneous data in a near-UV band is not con- 
sistent with the synchrotron component of the SEDs, while 
the (non-simultaneous) data in optical bands are well fit- 
ted by the synchrotron component. Since an one zone syn- 
chrotron model cannot satisfy both this near-UV point and 
a fiat spectrum in X-ray bands at the same time, this near- 
UV radiation is thought to be the contribution from an 
accretion disk. In ADAF models, the ICS component of 
high-temperature thermal electron s in the disk can mak e 
a spectral peak at th is band(e.g., iManmoto et al.l (|l997t ); 
iNemmen et al.l (|2010l )). The synchrotron component is not 
always dominated in all the bands below X-rays (see also 
Section U for possible contribution of the disk to X-rays). 

In order to check the uncertainty of Of, = 10 4 are 
considered in Fig[2] Since the cases of R = 10 14 and 10 15 
cm were ruled out, only SEDs for R = 10 16 and 10 17 cm are 
shown in the figure. Following the change of Of, the spec- 
tral edge of synchrotron radiation is shifted to lower energy 
by a factor of 10 2 . Also, the maximum energy of 7-rays de- 
creases, but this effect is smaller than the synchrotron edge 
because of the Klein-Nishina effect, which is the decrease of 
the IC cross-section compared with the Thomson scatter- 
ing. The detectability of the 7-rays by CTA for R = 10 16 is 
not largely changed. Note that both Figs[T] and [2] are con- 
sistent with the hard X-ray observations at present because 
the sensitivity of Swift BAT for 54 months is larger than 
the ex pected flux. Future obs ervations in hard X-rays, e.g., 
NeXT (|Takahashi et al] |2008). will constrain the maximum 
energy of electrons accelerated in the jets. 

Finally, we discuss SEDs in the case of S larger than 2.6. 




v[Hz] 

Figure 3. Same as Fig[T] but the cases of 5 = 10 and 6p = 10 4 . 



Fig[3] shows SEDs calculated on the assumptions of S = 10 
and Of = 10 4 . The other parameters are listed in Table [2] 
Since the high beaming factor boosts the flux of radiation, 
the energy density of radiation much lower than that in the 
case of Fig[T] is enough to reproduce the observed flux up 
to X-rays. Thus, the energy density of magnetic field in the 
blob is dominated compared to that of radiation even for 
R — 10 14 cm, but the cases of R ^ 10 15 ' 5 cm are not favored 
due to synchrotron self-absorption again. It is also domi- 
nated compared to the electron energy density, contrary to 
the cases of Figs [T] and [2] The smaller energy densities of ra- 
diation and electrons leads to much smaller ICS components 
than those in Fig[T] Thus, even long-term observations by 
Fermi LAT and CTA could not confirm these 7-rays. 



4 DISCUSSION & SUMMARY 

We have discussed possible 7-ray emission from jets in 
LLAGNs, and have also constrained several physical pa- 
rameters in our SSC model. The size of an emission region 
(blob) was limited to 10 16 < R < 10 17 5 cm by a monthly 
time scale variability and synchrotron self-absorption. This 
range is consistent with the variabilit y on typical time scale 
of a few days detected in LLAGNs |Anderson fc Ulvestadl 
l2005h . On the other hand, we have, so far, neglected a small 
(~ 10%) hours ti me scale (~ 1. 5 hr) v ariability in X-ray 
bands reported by lYounes et al.l (|2010l ). Here, we discuss 
the origin of this variability. 

This indicates the size of an emission region as 

R < S -^h±l = 5 x 10 14 5 3 AU Bil . 6h cm, (13) 
1 + z 

where At bs,i.sh = At b s /L5 hour is the time scale of 
the flux variability. This size is comparable with the 
Schwarzschild radius of the central black hole for S = 2.6 and 
is ~ 2 x 10 15 cm even for 5 = 10. These sizes are not favored 
due to the limitation of 7? if the X-rays are emitted from the 
blob. Thus, it is thought that X-ray radiation from an inner- 
most part of an accretion disk contributes to this variability. 
Theoretically, ADAF models can predict X-rays via thermal 
bremsstrahlung of ions in an inner hot part o f the disk (e.g., 
IManmoto et all (| 19971 ) INemmen et all j2010t )). This implies 
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that emission from the disk could contribute to the total X- 
rays by a small (but significant) fraction, even if a jet com- 
ponent is dominated from radio to X-ray bands suggested 
by the correlation between radio and X-ray luminosities. 

In this paper, a linear regime of synchrotron cooling of 
electrons, i.e., synchrotron cooling under constant magnetic 
field, was considered. Recently, non-linear radiative cooling 
of ele ctrons has been discussed |Schlickeiser fc Lerchdl2007l 
l200Sf ) . Non- linear models assume a constant equipartition 
parameter 77 _1 based on the success of spectral modeling 
of the observed blazars. Since the synchrotron and inverse 
compton cooling rates of electrons depend on the energy 
density of relativistic electrons, electron spectrum resultant 
from a solution of the kinetic equation of electrons is differ- 
ent from that derived from treatment in a linear regime even 
for a steady electron injection. Although whether equiparti- 
tion is realized in jets of LLAGNs is not clear observationally 
at present, it is an intriguing topic to investigate common 
physical features between blazars (strong radio galaxies) and 
LLAGNs. 

To summarize, we have discussed 7-ray emission from 
LLAGN jets in the framework of a SSC model on the as- 
sumption of radio and X-rays are dominantly produced from 
jets. The 7-rays are a direct probe of a jet component in ra- 
dio to X-ray bands without contamination from the other 
components, although the predicted flux is not large. Sev- 
eral observational results allowed us to constrain physical 
parameters in the emission region in jets. In the case of 
a beaming factor as low as that of parsec scale jets and 
R ~ 10 16 cm, CTA may detect the 7-rays in the near fu- 
ture and test the jet domination of radiation from LLAGNs. 
The determination of the respective contribution of disk and 
jet components will gives us a hint of a physical connection 
between a disk and relativistic jet in LLAGNs. 
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APPENDIX A: SYNCHROTRON SELF-COMPTON MODEL 

This appendix is dedicated to describe the detail of the synchrotron self-compton model used in this paper. We consider a 
steady state, that is a time-independent model. In addition, it is assumed that electrons and resultant synchrotron photons 
are distributed isotropically and magnetic field is randomly oriented in the blob. Below, the energies of electrons E e and 
photons Ej are measured in the unit of the electron mass, i.e., 7 = E e /m e c 2 and e 7 = E^/m e c 2 . 

The number density of photons in a blob frame, dn-Jde, is determined by a transport equation neglecting partial deriva- 
tives of time, 



1 

de 



(fik 



dtde 



"■7, ICS 
dtde 



(Al) 



t-7,csc 

where t 1>csc is the escape time scale of photons from a blob, d 2 n JtByn /dtde and d 2 n Ji ics/dtde are the production rate of 
synchrotron photons and ICS photons per unit volume from the electron distribution assumed in Eq[T] We assume an optically 
thin limit in this paper, that is t JyCBC = R/c. 

The production rate of synch rotron photons is calculated from the production rate of synchrotron photons from an 
electron in random magnetic field (|Crusius fc Schlickeiserll 198fj ), 



d 2 N ayn 3^3 a T cU B 2 
-Me; 7) = —r- 7. 9 



dtde 
where 



7r 7esem e c z 



Ki /3 {g)K 1/3 (g) 
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37 3 es ' 



es = 



7?n e c 



2-kUc 2r e U B 



(A2) 



(A3) 



and K p (x) is the modified Bessel function of p order. r e and h are the classical electron radius, Planck constant, respectively. 
Synchrotron radiation is accompanied by absorption in which a photon interacts wit h a charge in magnetic field. The absorption 
coefficient of synchrotron self-absorption is described in standard textbooks (e.g.. iRvbicki fe Lightmanl (| 19791 )) as 



8irce 



2nhc 



(h 



d 



d 2 N syll 2 _ 
dtde dj 



1 dn e 
7^~d7~ (7 ^ 



Thus, the production rate of synchrotron radiation is estimated as 



d n 



dtde 



^(e;7) + ^(e)«(e)c. 



The production rate of ICS is estimated as follows, 



d 2 mcs 



00 = 



d n. 



(0 



dtde dtde 
where #ics(7, e ) and Pics(e;7, e ') are 



J dy^-{-y)Rics(j, e) + c J de '-^( e ') J d 7 ^(7)- P 7,ics(e;7,e') 



flics (7, <0 = 



dfj, 



ctics(^,7, e ). 



(A4) 
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(A6) 



(A7) 



P7,ics(e;7. £ ') = 



dfj,- 



daics 



de 



(A8) 



H is the cosine of the angle between the momentums of incident photons and electr ons, cries is the cross-secti on of ICS in the 
case when an electron with the energy of 7 scatters a photon with the energy of e (|Coppi fc Blandfordlll99fj| ). 



CTICS (P, 7> <0 = 



(■, 4 8\. „ s 1 8 



1 



2(1 + : 



(A9) 



where x — 27e 7 (l — /3/x). daics/de is the differential cross-section of ICS in th e case when a photon with the energy of e' is 
upscattered into the energy of e by an electron with the energy of 7 |Leelll998h . 



daics , / v 3cr T 1 

— -, — ( £ ;7. £ >w = — 



+ 



+ 



(A10) 



Thus, flics ( 7 >e) an d Py,ics(e; 7, e') can be interpreted as angle-averaged cross-sections. 

After Eqs l A5I and IA6I are substituted, Eq |A2l is an integral equation cate gorized to the Fre dholm equation of the second 
kind. This equation can be numerically solved by the Nystrom method fe.g.. IPress et al] |l992j)). Note that dn-y/de includes 
higher order ICS photons, whose seed photons are ICS photons. Finally, photon flux observed at Earth is estimated as 



d 2 N^ 



dt bsde b: 



■(Cobs) = 8' 



V 



dn~, 



4-Kd L 2 t-, cac de 



[1 + z)e b s 



(All) 



where e Q b s = (1 + z) 1 8e and dh is the luminosity distance of the source. 



